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SECTION  1 
INTRODUCTION 


Parameters  of  interest  associated  with  the  passage  of  shock  wave 
over  a  dust  bed  include  the  mass  lofting  rate,  dust  velocity  profile, 
and  dust  mass  loading  profile*  The  data  on  these  and  similar  topics 
should  help  in  the  understanding  of  the  complicated  interaction  between 
the  gas  and  the  solid  phase  of  the  flow.  This  interaction  becomes  even 
more  complicated  in  the  presence  of  a  high  sound  speed  region. 

To  fully  justify  the  measurement  requirements,  a  comprehensive 
parametric  experimental  study  of  the  flows  was  required.  This  study  was 
beyond  the  scope  of  the  reporting  program.  However,  in  this  report,  the 
diagnostic  instrumentation  required  for  the  velocity  and  dust  density 
measurement  are  described  and  the  results  of  the  Initial  measurements 
are  reported. 

1.1  PREDICTION  OF  PARTICLE  TRAJECTORY. 

The  calculation  presented  here  was  performed  solely  to  obtain  a 
rough  estimate  of  the  particle  dynamics  in  order  to  guide  the  design  of 
the  experiment  and  to  help  analyze  the  measured  data.  To  determine  the 
measurement  location  of  the  multipoint  velocimeters ,  a  simple  numerical 
evaluation  of  the  particle  trajectory  was  made  by  solving  the  equations 
of  motion  for  an  isolated  particle.  The  primary  force  on  the  particle 
was  assumed  to  be  that  of  drag,  a  function  of  the  particle  Reynolds 
Number.  The  gravitational  force  was  also  included  in  the  equation  of 
motion  in  the  vertical  direction.  The  gas  velocity  profile  was  assumed 
to  be  known  and  taken  from  the  measurement  of  Reference  1.  Here,  data 
was  available  for  an  overpressure  of  23  psi  (158.6  k  Pa).  Initial 
trajectory  of  the  particle  into  the  boundary  layer  was  taken  from  the 
same  reference. 

The  particle  trajectory  calculation  was  first  carried  out  for  a 
2  mm  diameter  particle  for  which  the  data  were  available  (ref.  1).  The 
results  are  given  in  Figure  1.  It  shows  good  agreement  between  the 
experiment  and  the  computed  results. 
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Similar  calculations  for  overpressures  of  205  and  A 10  k  Pa  were 
carried  out.  In  the  absence  of  better  models,  the  gas  velocity  profiles 
were  assumed  to  scale  with  those  measured  for  the  overpressure  of  23  psi 
(168  k  Pa).  The  initial  oarticle  trajectory  angle  of  40  Degrees,  and 
Initial  particle  velocity  equal  to  2Z  of  the  freestream  gas  velocity 
were  assumed.  The  results  are  shown  in  Figures  2  and  3.  Figure  2  shows 
the  trajectory  of  the  particles  during  the  initial  phase  of  the  lift 
off.  The  smaller  particles  lose  their  initial  momentum  rapidly  and 
follow  the  gas  stream.  The  larger  particles,  due  to  their  larger 
initial  momentum,  penetrate  the  gas  boundary  layer.  The  results  show  the 
classic  mechanism  of  aerodynamic  sorting  of  the  particles. 

It  should  be  pointed  out  that  the  smaller  particles  are  more 
influenced  by  the  turbulence  within  the  boundary  layer  (not  modeled 
here),  and  hence  are  apt  to  be  more  uniformly  dispersed  within  the  flow 
field.  Hence,  a  uniform  sorting  of  the  particles  according  to  their 
size,  as  predicted  here,  should  not  be  expected.  On  the  other  hand,  the 
larger  particles  are  less  affected  by  the  velocity  fluctuations  so  their 
distributions  are  probably  more  accurately  presented  here. 

Figure  3  shows  the  trajectories  of  different  particles  and  their 
expected  elevation  as  a  function  of  time.  The  position  of  the 
measurement  stations  are  also  included.  Based  on  this  model,  the  length 
of  the  dust  bed  must  be  at  least  17  meters  so  as  to  expect  a  100  micron 
particle  to  reach  a  height  of  20  centimeters,  our  highest  measurement 
locations  at  CERF  6'  shock  tube.  These  numbers  were  used  as  a  guideline 
for  the  CERF  6'  shock  tube  measurements. 

An  estimate  of  the  particle  relaxation  tiue  for  two  overpressures 
of  205  and  410  psi  were  made  and  shown  in  Figure  4.  It  shows  that 
first,  the  overpressure  does  not  greatly  influence  the  relaxation  period 
of  the  particles,  and  secondly,  a  relaxation  time  of  the  order  of  10  ms 
(for  predicted  velocity  jump  across  the  shock  wave)  is  predicted  for  100 
micron  particles.  This  information,  used  as  an  order  of  magnitiude 
information,  helps  determine  the  minimum  measurement  time  required  for 
any  given  particle  size  to  reach  equilibrium  within  the  carrying  flow. 


SECTION  2 
INSTRUMENTATION 


The  components  of  the  experimental  subject  matter  are  briefly 
described  herein.  These  Include,  the  principle  of  operation  for  the 
measurement  technique,  and  the  description  of  the  apparatus  designed  for 
the  measurments.  In  each  case,  the  shock  tube  used  in  the  experiment  is 
briefly  described,  and  then  the  apparatus  designed  for  the  measurement 
is  reviewed. 

2.1  PRINCIPLE  OF  OPERATION. 

Direct  measurement  of  the  particle  velocity  was  obtained  using 
Laser  Doppler  Velocimetry  (LDV)  [Ref.  2  and  3].  The  apparatus  generally 
consists  of  a  coherent  source  of  light  (laser)  included  in  a  trans¬ 
mitting  optics,  and  receiving  optics  as  depicted  in  Figure  5.  Two 
coherent  beams  are  focused  at  the  measurement  point,  known  as  the  probe 
volume,  resulting  in  the  formation  of  the  interference  fringes.  The 
distance  between  fringes  (5)  is  determined  by  the  beam  angle  (9)  and  the 
laser  beam  wavelength  (X)  as  given  by 

0  -  X/2  sin( 9/2)  (1) 

When  a  particle  with  a  velocity  U  crosses  the  probe  volume,  it  scatters 
light  proportional  to  the  local  intensity  of  the  light.  The  receiving 
optics  collects  portion  of  the  light  and  transforms  it  into  an  electri¬ 
cal  signal,  known  as  Doppler  burst.  The  period  of  the  burst  (t)  is  a 
function  of  the  fringe  spacing  and  the  instantaneous  velocity  of  the 
particle,  U.  Through  direct  digitization  of  the  signal  and  post  pro¬ 
cessing,  the  particle  velocity  may  be  obtained  from  the  relation 

U  -  ( <5/ t)  (2) 

Simultaneous  measurement  of  the  opacity  of  the  measurement  media  is  made 
with  the  use  of  a  photodiode  as  a  beam  stop.  The  reduction  of  the  total 
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optical  energy  collected  on  the  photodiode  can  theoretically  be  related 
to  the  increase  in  the  dust  loading  in  the  shock  tube  test  section. 
Here,  assumptions  are  made  that  the  scattering  characteristics  of  Che 
dust  particles  are  uniform;  that  the  amount  of  light  scattered  is 
proportional  to  the  scattering  cross  section  of  the  particles;  and  that 
the  effect  of  the  window  contaminations  are  negligible. 


2.2 


TRW  4  X  4  INCH  SHOCK  TUBE  . 


A  schematic  of  the  shock  tube  at  TRW  is  shown  in  Figure  6.  It 
has  a  cross  section  of  3.5  by  3.5  Inches  (8.9  by  8.9  cm)  with  the  floor 
replaced  by  a  dust  bed  nominally  15  cm  inches  deep.  The  dust  particle 
size  range  was  1  to  400  microns,  with  a  mass  mean  of  120  microns.  The 
shock  tube  configuration  characteristics  are  summarized  in  Table  1. 
Optical  access  to  the  flow  field  was  made  through  two  15  cm  diameter 
windows.  The  test  duration  was  5  ms,  limited  by  the  dust  bed  leading 
characteristic  (Ref.  4,  5].  The  shock  overpressure  at  the  test  section 
was  nominally  at  210  k  Pa. 


Laser  Velocimeter.  Measurement  of  the  dust  particle  velocity  was  made 
using  a  single  component,  fringe  mode.  Laser  Doppler  Velocimeter  as  shown  in 
Figure  7.  A  4-watt  Argon  Ion  laser  operating  in  single  frequncy  mode 
was  used.  The  light  beam  was  focused  onto  a  bragg  cell  used  both  as  a 
beamsplitter  and  for  frequency  shifting.  The  beams  were  further 
collimated  through  a  second  lens,  and  crossed  over  at  the  probe  volume 
by  the  transmitting  lens.  The  scattered  light  was  collected  in  back 
scatter  mode  onto  the  end  of  a  fiber  optic  cable  which  carried  the 
collected  scattered  light  to  a  suitable  photomultiplier.  A  photograph 
of  the  optical  system,  Including  the  shock  tube  is  shown  in  Figure  8. 

An  important  criteria  to  the  success  of  the  measurement  was  to 
maximize  the  number  of  samples  collected  during  each  time  window.  Since 
the  number  of  scattering  centers  could  not  be  increased,  as  they  were 
provided  by  the  dust  bed,  the  data  collection  was  designed  to  maximize 


Table  1 


TRW  4“  Shock  Tube  Conditions 


Nominal  Test  Overpressure 

Nominal  Test  Time 

Test  Section  Dimension 

Test  Section  Length 

Dust  Bed  Material 
40-48%  Porosity 

Specific  Gravity 
1-5  (noncompacted) 

1-74  (compacted) 

Particle  Size 

average  ~  120  urn 


210  k  Pa 
5  ms 

8.9  cm  x  8.9  cm 
1.57  m 

White  Sands  Natural 
2-9  (solid  only) 

d  <  420  pm,  mass  mean 


the  acceptance  rate  of  Doppler  bursts.  The  pedestal  of  the  amplified 
signal  was  first  removed  through  a  high  pass  filter.  The  signal  was 
then  mixed  down  with  a  30  MHz  signal  to  enable  the  signal  to  be  digi¬ 
tized  at  a  rate  of  100  MHz.  With  the  available  total  memory  of  512K  a 
rerl  time  recording  duration  of  5.12  ms  was  achieved.  The  recorded  data 
were  transferred  to  recording  disks  for  archiving  and  retrieval.  The 
block  diagram  of  the  data  acquisition  system  is  shown  in  Figure  9. 

2.3  CERF  6'  SHOCK  TUBE. 

The  6'  Shock  Tube  located  at  Civil  Engineering  Research  Facility 
(CERF)  is  made  of  6'  (1.97  m)  diameter  steel  sections  shown  in 
Figure  10.  The  velocity  measurements  were  made  at  Section  G,  71  m  from 
the  opening  in  the  driver  side  and  approximatly  23  m  into  the  dusty 
section. 

The  instrumentation  was  inside  a  concrete  bunker  placed  under¬ 
neath  the  shock  tube  as  shown  in  Figure  11.  The  optical  probe  was 
placed  over  a  traversing  mechanism,  Figure  12,  and  raised  into  the  shock 
tube  through  an  opening  at  the  bottom  of  the  tube.  The  probe  inside  the 


shock  tube  was  protected  by  a  blast  wing  shown  in  Figures  13  and  14. 

The  window  plate  was  designed  such  that  the  measurement  height,  defined 
by  the  position  of  the  optical  probe,  could  be  continuously  varied. 
Figure  13  shows  the  photograph  of  the  concrete  bunker  used  in  this 
experiment.  The  shock  tube  configuration  is  given  in  Table  2. 

Table  2.  Cerf  6'  Shock  Tube  Conditions. 

Nominal  Test  Overpressure  210  k  Pa 

Nominal  Test  Time  5  ms 

Test  Section  Length  29.5  m 

Test  Bed  Material  McCormick  Ranch  Natural 

Particle  Size  d  <  5000  y 

mean  100  y 

2.3.1  Multipoint  Laser  Velocioeter  . 

A  four-point  simultaneous  velocimeter  was  designed  for  this  test 
using  a  4-watt  Argon  Ion  laser.  Again,  a  bragg  cell  was  employed  as  a 
beamsplitter  and  for  frequency  shift,  as  shown  in  Figure  16.  The  two 
coherent  beams  of  equal  Intensity  were  divided  into  four  separate  LDV's, 
through  a  series  of  beamsplitters.  Approximately  25%  of  the  laser  power 
was  directed  to  each  station.  The  transmitting  optics  consisted  of 
turning  mirrors  and  focusing  lenses.  The  receiving  optics  were  also 
placed  on  the  same  support  beam  as  the  transmitting  optics,  and  focused 
the  scattered  light  onto  four  fiber  optics  cables.  The  details  of  the 
optical  head  assembly  is  shown  in  Figure  17.  The  optical  probe  and  a 
close  up  of  the  optical  head  assembly  are  shown  in  Figures  18  and  19. 

2.3.2  Data  Acquisition. 

The  system  diagram  for  the  four  component  laser  velocimeter  is 
shown  in  Figure  20.  The  four  fiber  cables  transmit  the  collected  light 
to  Photomultiplier  Tubes.  The  signals  are  then  amplified  and  run 
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through  a  multiplexer  which  transmits  one  signal  according  to  a 
preprogramed  schedule.  This  multiplexing  is  done  in  order  to  conserve 
the  digital  memory.  The  output  of  the  multiplexer  is  then  mixed  down 
and  recorded  using  a  100  MHz  transient  recorder.  The  entire  data 
acquisition  and  trigger  and  firing  sequence  is  controlled  by  an  IBM 
PC.  A  photograph  of  the  electronic  instrumentation  is  shown  in 
Figure  21. 

2.3.3  Data  Analysis. 

An  important  aspect  of  the  data  acquisition  of  a  laser 
veloclmeter  system  is  the  method  by  which  the  signal  is  processed. 

Here,  the  velocity  information  (signal)  should  be  separated  from  the 
anomalies  (noise)  in  some  form  or  fashion.  In  the  present  system,  since 
the  unedited  data  were  recorded,  various  analysis  algorithms  were 
applied  to  an  identical  set  of  data  and  the  relative  effectiveness  of 
the  techniques  were  evaluated.  Two  of  the  major  techniques  which  were 
implemented  are  briefly  discussed  here. 

2.3.4  Three  Level  Comparator  and  a  Phase  Noise  Detector. 

A  software  version  of  the  "three  level  comparator  technique"  was 
Implemented.  This  technique  has  been  successfully  used  in  commercially 
available  analogue  counters.  Here,  positive  and  negative  voltage 
thresholds  are  chosen.  Each  filtered  burst  is  evaluated  by  the  sequence 
at  which  the  signal  crosses  the  two  thresholds  and  the  zero  voltage 
line.  Here,  the  criteria  is  that  at  least  one  digitized  point  lies 
between  two  successive  zero  crossing  points.  In  practice,  the  threshold 
values  are  chosen  to  be  above  the  background  noise,  and  yet  to  be 
smaller  than  the  signal  voltages.  The  effectiveness  of  the  technique 
depends  on  the  choice  of  the  threshold.  In  our  case,  as  we  had  access 
to  the  actual  signal  which  is  going  to  be  processed,  the  level  of  the 
background  noise  during  each  time  window  was  evaluated  and  the  values  of 
the  threshold  were  automatically  optimized.  Figure  22  shows  a  typical 
burst  processed  using  this  technique. 
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A  criteria  for  detection  of  phase  noise  in  the  signal  was 
developed.  Phase  noise  usually  is  generated  when  more  than  one  particle 
is  present  at  the  probe  volume.  Here,  the  time  rate  of  change  of  the 
signal  was  used  as  the  test  vehicle.  In  the  absence  of  phase  noise,  the 
first  derivative  of  the  signal  between  any  two  successive  zero  crossing 
changes  sign  no  more  than  two  times.  A  signal  with  more  than  two 
changes  of  the  sign  of  its  derivative  between  successive  zero  crossings 
was  then  rejected. 

An  additional  criteria  comparing  the  calculated  period  of  two 
different  numbers  of  cycles  of  the  s use  burst  was  also  Included  in  the 
discrimination  software.  This  provided  an  additional  test  for  the 
accuracy  of  the  discriminating  technique.  Through  the  software,  the 
effects  of  each  discriminating  technique  were  evaluated.  It  was 
determined  that  the  combination  of  the  above  mentioned  techniques,  with 
as  few  as  three  full  cycles  of  the  burst,  was  sufficient  for  processing 
the  signal  obtained  in  the  shock  tube. 

2.3.5  Burst  Correlation/Prequency  Spectrum  Technique. 

Another  technique  which  was  found  to  be  more  powerful  for  the 
cases  of  low  slgnal-to-noise  ratios  and  limited  sample  Is  described 
here.  The  particle  velocity  is  assumed  to  be  constant  In  every  time 
window  of  128  microseconds.  Within  each  time  window,  bursts  are 
recognized  by  a  voltage  threshold  test.  The  acceptable  bursts  are 
recognized  through  their  respective  power  spectra.  Fast  Fourier 
Transformation  is  used  and  the  signal  power  distribution  in  the 
frequency  domain  is  obtained.  A  criteria  is  used  here  to  distinguish 
between  acceptable  and  unacceptable  signals.  The  acceptable  bursts  are 
further  smoothed  through  autocorrelation.  Finally,  the  mean  velocity 
and  standard  deviation  are  obtained  by  Fourier  Transformation  of  the 
mean  autocorrelation  of  the  bursts  during  of  time  window.  Figure  23 
shows  a  typical  burst  processed  using  spectrum  analysis. 

Figure  24  exhibits  three  typical  "bursts”  reproduced  from  the 
shock  tube  data  set.  They  represent  a  strong  signal,  a  marginal  signal, 


and  what  appears  to  be  a  very  weak  signal  or  possibly  a  random  noise. 
When  a  "three  level  comparator"  discrimination  is  applied  to  the  above 
traces,  only  the  first  trace  Is  accepted  and  the  remainder  are 
rejected.  On  the  other  hand,  when  the  power  spectrum  of  the  signals  are 
calculated,  as  shown  in  Figure  24,  the  frequency  content  and  the 
strength  of  the  signal  for  each  burst  is  evaluated.  In  this  case,  the 
dominant  frequency  in  all  three  bursts  are  close  to  each  other,  while 
the  signal  strength  is  different. 


SECTION  3 

MEASUREMENT  MATRIX 


Limited  number  of  shock  tube  entries  were  made  during  the  course 
of  the  present  contract.  The  measurements  were  conducted  in  dusty  non- 
precursed  flows.  The  run  conditions  for  the  two  shock  tubes  vrere 
generally  different  and  no  attempt  is  made  here  to  correlate  data  beyond 
the  overall  features  of  the  dusty  boundary  layer  flows.  More 
comprehensive  Investigation  of  the  flows  is  curently  being  carried  out 
using  the  instrumentation  described  in  this  report  through  our  present 
contract  (DNA001-84-C-0116)  with  the  Defense  Nuclear  Agency  (DNA). 

3.1  TRW  4  X  4  INCH  SHOCK  TUBE  . 

Measurements  were  obtained  at  a  test  section  located  3.6  m  from 
the  leading  edge  of  the  dust  bed.  The  laser  veloclmeter  was  placed  on 
an  isolation  table.  The  height  of  the  probe  volume  was  adjusted  through 
elevation  of  the  entire  table.  The  shock  tube  operated  by  pressurizing 
the  driver  until  bursting  an  aluminum  diaphragm.  The  general  flow 
parameters  of  the  flow  field  are  given  in  Table  1.  The  dust  bed 
preparation  consisted  of  running  the  tube  a  few  times  to  achieve 
complete  compactness  of  the  dust  bed.  The  surface  of  the  dust  bed  was 
smoothed  prior  to  each  shot.  The  flow  conditions  were  generally  kept 
uniform  for  subsequent  shots. 

The  shock  tube  driver  was  operated  for  a  nominal  210  k  Pa 
overpressure  at  the  test  section  (with  no  dust).  The  strength  of  the 
shockwave  however,  was  reduced  as  it  passed  over  the  dust  bed.  The 
pressure  time  history  at  the  measurement  location  is  shown  in  Figure  25. 
It  shows  a  square  wave  profile  for  a  duration  of  approximately  5  ms. 

The  limiting  wave  was  the  compression  from  the  leading  edge  of  the  soil 
bed.  The  peak  over  pressure  was  measured  to  be  roughly  138  k  Pa. 


3.2  CERF  6’  SHOCK  TUBE  . 

Measurements  were  obtained  at  a  test  section  located  70  meters 
from  the  start  of  the  dust  bed.  The  height  of  the  probe  was  varied 
using  the  traversing  mechanism  (Figure  12).  The  general  flow  parameters 
are  given  in  Table  2.  The  dust  bed  was  compacted  by  a  roller  before 
each  shot,  and  the  driver  was  kept  in  the  same  configuration  with  the 
same  charge  level  to  maintain  the  same  flow  conditions  for  each  shot. 

The  shock  tube  driver  produced  a  nominal  overpressure  of  210  k  Pa 
at  the  measurement  position.  The  time  history  at  this  position  is  shown 
in  Figure  26.  It  shows  that  the  pressure  declines  approximately  35  k  Pa 
for  the  5  ms  duration. 


row 


SECTION  4 
RESULTS 


Presentation  of  the  results  of  the  measurement  are  made  in  this 
chapter.  The  data  are  grouped  in  three  categories: 

1.  Dust  velocity  data  at  TRW  4"  shock  tube. 

2.  Dust  concentration  data  at  TRW  4"  shock  tube. 

3.  Dust  velocity  data  at  CERF  6*  shock  tube. 

In  each  case,  the  data  acquisition  was  initiated  by  the  passage  of  the 
shock  wave  at  the  measurement  point  and  stopped  after  the  memory  bank 
was  filled.  No  data  compression  was  attempted.  With  512k  bytes  of 
memory  available,  the  total  measurement  duration  in  each  case  was 
limited  to  about  5  ms,  coinciding  roughly  with  the  run  duration  of  the 
shock  tubes. 

4.1  DUST  VELOCITY  DATA  AT  TRW  4"  SHOCK  TUBE  . 

During  each  run  of  the  shock  tube,  only  one  component  of  the 
particle  velocity  was  obtained.  The  dust  surface  preparation  and  the 
driver  conditions  were  repeated  for  subsequent  runs.  The  shock  tube 
conditions  and  the  static  overpressure  profiles  are  given  in  Table  l. 

Figures  27  to  30  give  the  velocity  data  for  the  measurement 
heights  of  6  mm  (0.25  inch)  to  25.4  mm  (1  inch)  above  the  initial  dust 
level.  In  each  case,  the  data  obtained  at  the  same  height  for  different 
runs  are  presented  together.  The  velocity  data  at  6  mm  (.25  inch)  is 
shown  in  Figure  27.  The  first  arrival  of  particles  at  the  measurement 
location  was  about  .25  ms  after  the  passage  of  the  shock  wave.  This 
arrival  time  was  consistent  for  all  runs.  During  the  subsequent  1.5  ms, 
an  Increase  in  the  particle  velocity  is  observed.  This  represents  the 
response  of  the  particles  to  the  step  function  in  the  gas  velocity 
behind  the  moving  shock  wave. 


The  vertical  component  of  the  dust  particle  velocity  at  the  same 
measurement  height  is  also  given  in  Figure  28.  Measurement  consistency 
of  the  runs  presented  are  very  good.  The  normal  velocity  component 
shows  a  slow  change  in  the  angularity  of  the  flow  of  the  dust 
particles.  At  1.5  ms  after  passage  of  the  shock,  the  flow  angle  is 
measured  to  be  about  10  degrees.  This  angle  indicates  the  development 
of  a  boundary  layer  near  the  dust  surface. 

The  development  of  the  boundary  layer  is  more  apparent  from  the 
velocity  data  recorded  at  the  height  of  6  mm  (same  as  Figure  27),  but 
recorded  for  4  ms  as  presented  in  Figure  31.  The  dust  velocity  shows  a 
gradual  decline  after  1.5  ms,  representing  the  growth  of  the  boundary 
layer  beyond  the  measurement  height. 

Similar  data  are  obtained  at  heights  of  13  mm  (.5inch),  19  mm 
(.75inch),  and  25  mm  (1.0  inch),  and  are  shown  in  Figures  28  to  30.  The 
overall  behavior  of  the  dust  particles  are  similar.  This  similarity  is 
further  demonstrated  in  Figure  32,  where  all  the  measurements  are  drawn 
together.  It  clearly  shows  the  initial  particle  relaxation  period  of 
approximately  2  ms.  The  early  time  ( t  <  1  ms)  velocity  development 
behind  the  shockwave  show  a  plug  flow  type  of  velocity  field.  The 
scatter  of  the  data  in  this  region  may  be  attributed  to  the  variation  in 
the  particle  size  distribution,  as  well  as  variation  in  the  shock  tube 
condition.  The  development  of  the  velocity  gradient  (and  hence,  the 
velocity  boundary  layer)  is  also  apparent.  Based  on  the  data  of 
Figure  32,  it  is  estimated  that  the  velocity  boundary  layer  reaches  the 
highest  measurement  height  of  25  mm  (1  inch)  at  about  4  ms. 

From  the  velocity  time  histories  of  proceeding  figures,  two 
boundary  layer  type  profiles  may  be  defined.  One  is  determined  by  the 
arrival  time  of  the  dust  particle  at  each  of  the  measurement  heights  as 
shown  in  Figure  33.  The  second  profile,  also  shown  in  the  same  figure, 
is  defined  by  the  points  in  Figure  32,  for  which  the  time  rate  of  change 
of  velocity  becomes  negative.  This  suggests  that  in  dusty  flows,  the 
customery  definition  of  the  boundary  layer  must  be  augmented  to 
distinguish  between  it  and  the  boundary  layer  defined  by  the  presence  of 
the  dust  particles. 


An  important  piece  of  data  used  for  the  evaluation  of  the 
computation  prediction  of  the  dusty  flow  is  the  velocity  profile.  Here, 
the  velocity  profile  obtained  at  4  ms  behind  the  shock  wave  is  plotted 
in  Figure  34.  The  height  of  the  "velocity"  boundary  layer  was  assumed 
to  be  25  mm  (1  inch).  It  was  apparent  that  four  points  were  not 
sufficient  to  clearly  define  the  velocity  profile.  But  it  was  also 
recognized  that  within  the  constraints  of  the  measurement  environment 
and  the  uncertainty  of  the  dust  bed,  measurement  closer  to  the  surface 
was  not  possible.  Figure  34  also  shows  the  velocity  data  to  best  fit  a 
1/4  or  1/5  power  law  profile. 

4.2  DUST  CONCENTRATION  DATA  AT  TRW  4"  SHOCK  TUBE  . 

As  shown  in  Figure  5  of  Section  2,  a  photodiode  was  used  as  the 
laser  light  beam  stop  when  measurement  at  TRW  4"  shock  tube  was  made. 

The  output  of  the  photodiode  registered  the  total  laser  power 
transmitted  through  the  shock  tube  test  section.  Figures  35  to  38  show 
the  raw  data  obtained  at  the  measurement  heights.  The  reduction  in  the 
transmitted  laser  power  after  the  passage  of  the  shock  wave  is  caused  by 
the  presence  of  the  scattering  centers  (dust  particles)  along  the  path 
of  the  beams  in  the  shock  tube.  This  response  Includes  the  buildup  of 
the  dust  on  the  window  which  at  its  highest  is  estimated  to  scatter 
about  20 %  of  the  incident  laser  power.  It  was,  however,  decided  that 
due  to  the  uncertainty  of  the  rate  of  the  dust  buildup  on  the  window 
during  the  first  few  milliseconds,  no  correction  will  be  applied  to  the 
measured  data.  The  combined  plots  of  the  extinction  curves  are  given  in 
Figure  39. 

Data  on  the  particle  size  distribution  at  the  measurement 
location  is  required  to  reduce  the  laser  light  extinction  to  dust 
concentration  data.  This  information  was,  however,  not  available  for 
all  time  as  independent  measurement  of  the  particle  size  was  not 
measured.  (This  measurement  is  presently  being  carried  out  at  the  TRW 
17"  shock  tube).  There  were,  however,  regions  in  the  early  measurement 
time  history  where  it  was  determined  that  the  particle  size  distribution 
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to  be  roughly  between  1  and  10  micrometers.  This  time  region  was 
bounded  by  the  onset  of  the  laser  extinction  light,  and  the  arrival  of 
the  larger  than  10  micron  particles  determined  from  the  velocity 
measurements  and  shown  in  Figure  40.  Based  on  the  uniform  weight 
distribution  (cube  root  particle  size/number  density  distribution),  the 
average  particle  size  for  this  time  period  was  4  micrometers. 

Based  on  the  assumed  particle  size  distribution,  the  theoretical 
calibration  curve  of  Figure  41  was  used  to  convert  laser  light 
extinction  to  dust  concentration.  The  results  are  given  in  Figure  42. 
In  the  absence  of  independent  data  on  particle  size  distribution  beyond 
1  ms,  the  determination  of  the  dust  density  beyond  this  time  was  not 
attempted. 

4.3  DUST  VELOCITY  DATA  AT  CERF  6’  SHOCK  TUBE  . 
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Fewer  measurements  were  attempted  at  the  6'  shock  tube  during  the 
reporting  contract.  The  results  of  the  measurements  for  a  series  of 
nonprecursed  dusty  flows  are  reported  in  Figures  43  to  46.  Here,  the 
velocity  measurements  represent  the  same  trend  as  those  obtained  in  the 
smaller  shock  tube. 

Figure  47  shows  the  velocity  data  for  all  the  measurement 
heights.  The  development  of  a  velocity  boundary  layer  is  apparent  as 
early  as  2  ms  after  the  passage  of  the  shock  front.  Figures  48  and  49 
present  the  velocity  profiles  obtained  from  the  velocity  data.  Also 
from  the  velocity  data,  Figure  50  shows  the  time  of  arrival  of  the  dust 
particles  at  the  various  measurement  heights.  No  measurement  of  the 
vertical  velocity  components  were  made  in  this  facility. 

Preliminary  measurements  of  the  dust  velocity  in  precursed  dusty 
flow  was  made.  The  dust  velocities  were  found  to  be  generally  higher 
than  those  obtained  for  the  nonprecursed  flows. 


SECTION  5 
CONCLUSIONS 


The  measurements  at  both  TRW's  4x4  inch  shock  tube  and  the  CERF 
6*  shock  tube  were  obtained  at  various  points  above  the  dust  bed.  It 
was  shown  that  the  particle  velocity  just  behind  the  shock  wave  was  not 
a  strong  function  of  the  measurement  height.  The  results,  however,  show 
the  development  of  a  boundary  layer  for  both  test  sites  at  a  laser  time 
behind  the  shock  wave.  For  the  TRW  data,  the  velocity  profile  was 
approximately  similar  to  that  of  a  turbulent  boundary  layer  with  a  rough 
surface,  following  a  power  law  with  exponent  values  of  1/4  or  1/5.  The 
results  at  the  CERF  6'  shock  tube  were  similar  with  an  exponent  value 
approximately  1/4. 

An  attempt  was  made  to  estimate  the  density  of  the  particles  at 
some  positions  in  the  flow  for  the  4  inch  shock  tube  at  TRW.  Simulta¬ 
neous  measurements  of  the  laser  light  extinction  was  made.  This  infor¬ 
mation  along  with  assumptions  on  the  particle  size  distribution 
generated  the  data  on  the  particle  concentration  behind  the  moving  shock 
wave. 
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Figure  9.  Block  Diagram  of  the  Data  Acquisition  System 


Figure  11.  Laser  Velocimeter  Installation  at  CERF  6'  Shock  Tube 
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Figure  20.  System  Diagram  for  Laser  Velocioeter  at  CERF  6'  Shock  Tube. 
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Figure  31.  SCreamwise  Particle  Velocity  at  TRW  4"  Tube, 
Shock  Wave,  Y  ■  6  an. 
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Figure  35.  User  Light  Extinction,  Y  -  6 
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Figure  39.  laser  light  Extinction  Across  the  Shock  Tube  for  All 
Heights* 
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Figure  40.  Dust  Particle  Size  Group 
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Figure  43.  Particle  Velocity  Tim  History  et  CERF  6’  Shock  Tube, 
Y  -  12  aa. 
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Figure  48.  Velocity  Profile  it  2  ■  Behind  Shock  Have 
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Figure  49.  Velocity  Profile  et  4  wm  Behind  Shock  Wave 
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